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The distance to IC4499 from 
K-ban6 photometry of 32 RR Lyrae stars* 
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Abstract. Single epoch A'-band photometry of 32 RR Lyrae stars in the young galactic globular cluster IC4499 is presented. 
The mean magnitudes have been estimated using the AT-band template light-curves of Jones, Carney and Fulbright 1 1996 1. We 
derive an independent estimate of the distance and reddening for the cluster based on the new (Mk) -log P -[Fe/H] relation from 
Bono et al.(2003 ) which has a zero-point in agreement with the canonical LMC distance of (m - M)o = 18.5. For an assumed 
metallicity of [Fe/H] = -1.65 ± 0.1 we find (m - M)q = 16.47 ± 0.04 ± 0.06 (random and systematic errors respectively) and 
E(B - V) = 0.24 ± 0.03 in good agreement with the findings of Walker and Nemec <1996> based on optical data. The distance 
estimate is very robust to uncertainties in the reddening estimate but naturally depends directly on the adopted zero-point of the 
(Mk) -log P -[Fe/H] relation. 
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1. Introduction 

The globular cluster IC4499 (RA: 15''00'"19^ DEC: 
-82°12'50"; 2000.0) is among the youngest galactic globular 
clusters (Ferraro et al. 119951 Fusi Pecci et al., [p95 1. It is 
of OosterhofF type I and is thus likely to have zero or even 
retrograde rotation around the Galactic center (van den Bergh 
119931 Lee and Carney 1999)- It is a fairly low density cluster 
with one of the highest numbers of RR Lyrae stars per unit 
luminosity (Suntzeff, Kinman and Kraft I199I> containing 
about 100 RR Lyrae stars. These extreme characteristics 
suggests that the cluster might be a late addition to the Galaxy, 
which makes the cluster particularly interesting to study. 

Optical color-magnitude diagrams have been presented by 
Sarajedini (I1993> and Ferraro et al. iI995> . the latter extend- 
ing two magnitudes below the main-sequence turn-off. The 
RR Lyrae stars have been the subject of several investiga- 
tions (Cout ts et a l. [T9751 Clements et al. [T979I [T986I Walker 
and Nemec 1996^(WN96 in the following)) since their discov- 
ery (Fourcade and Laborde, 1969 and Fourcade et al., 1974 1. 
These investigations have revealed that the RR Lyrae stars span 
the full range of pulsation modes cuiTently known, including 
stars exhibiting double mode pulsations, which can be used to 
determine the mass of these stars. WN96 has even suggested 
that some of the stars might pulsate in the second overtone. 
The range of the stellar masses has been shown by Sarajedini 



JI993> to be very limited as judged from the double mode pul- 
sators. 
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The reddening towards IC4499 is fairly high, E(B - V) - 
0.22 + 0.02 (WN96) and thus a distance estimate which is in- 
sensitive to reddening will provide an important additional con- 
straint on the distance to this cluster. 

The (Mk) - log P relation for RR Lyrae stars, pioneered by 
Longmore et al. J 1 9861 1 and refined further by Longmore et al. 
( 1990), is a very powerful tool for distance determination. The 
method is theoretically based on the pulsation equation from 
van Albada and Baker J197II I which links the period to the lu- 
minosity of an object. Bono et al. ( I2001» . and most recently 
Bono et al. ( 2003 ) (B03 hereinafter) has refined the method us- 
ing pulsational stellar models to include a small but significant 
metallicity term, and we will use this calibration in the follow- 
ing. It has been shown by Butler ( 2003 1, through a comparison 
of inner and outer region RR Lyrae stars in M3, that the method 
is insensitive to environmental effects. It is also rather insensi- 
tive to errors in both [Fe/H] and reddening, and as the ampli- 
tudes of RR Lyrae stars in the ^T-band are smaller than about 
0.3 mag it is feasible to apply the relation to single epoch data 
if a sufficiently large number of stars can be observed. Using 
recent ephemerides for the IC4499 stars we can even phase the 
observations and use template light-curves to improve our esti- 
mate of the mean /T-magnitudes. 

2. The Data 

We used the IRAC-2B near-IR camera (IVloorwood et aL ll992> 
mounted at the ESO-MPI 2.2m telescope at La Silla on the 
night Sep. 17/18, 1995. The detector was a 256 by 256 pixel 
NICMOS-3 array which was read out using the so called "non- 
destructive read" mode. The lens LC, which provides a pro- 
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jected pixel size of 0.5 arcsec was employed, and the seeing 
was typically 1.3 arcsec FWHM. 

Two adjacent fields (labeled A and B in the following) in- 
cluding the center of IC4499 were selected on the basis of 
the charts in Clement et al. J1986> to cover a large number of 
known RR Lyrae variables as well as a good range in period. 

The observations were performed in a regular pattern of 
two dithered target frames interspersed by two dithered sky 
frames (several arc-minutes away) followed by the same se- 
quence again, after offsetting the telescope by a few arc- 
seconds. Each frame was made of six 10 seconds integrations 
which were averaged on-line by the pre-processing system. For 
field A, 16 individual frames were obtained through the K' fil- 
ter and for field B, 11 such frames were obtained giving a to- 
tal of 16 and 11 minutes, respectively, of on-target integration 
time. The mean HJD for the two datasets is 2449978.5198 and 
2449978.5661 respectively. 

3. Data reduction 

3.1. Pre-processing 

All the data were corrected for non-linearity using the 
ctio.irlincor command within IRAF'. Polynomial coef- 
ficients of coeffl^l.O, coefr2= 3.808477 x 10"^ coefr3= 
5.37237 X 10"^ were adopted on the basis of the reported non- 
linearity for the array (Lidman. ll995> . 

Flatfields were produced as the difference between expo- 
sures of the white screen in the dome with the quartz lights 
on and off, respectively. The IRAC-2 camera is known to suf- 
fer from parasitic light in the flatfields at the 5% level, and a 
standard star was mapped in a 4 by 4 grid across the array on 
the night of the IC4499 observations as well as on the previ- 
ous night in order to determine the necessary correction for 
the flat fields. The counts in each position were fitted by a two 
dimensional second order polynomium for each night. The re- 
sults were in good agreement, to within 3%, and the averaged 
correction was applied to the flat field frame. 

Sky frames were computed as the average of four dithered 
sky exposures obtained before and after the science frame. The 
"avsigclip" algorithm of the IRAF imcombine command with 
a threshold value of cr = 5 was employed to eliminate stars 
in the sky fields. For each target frame the corresponding sky 
frame was subtracted and bad pixels were removed by linear in- 
terpolation. The resulting targets frames were finally averaged 
to produce the final target frames. 

3.2. Standard star photometry 

Fourteen standard star measurements were performed on the 
night Oct. 1/2, 1995, to determine the color terms of the cam- 
era in the J and K' filters and the airmass coefficients. Each 
star was observed 4 times with the star imaged in four different 
positions on the array. The stars were selected from the list of 
Carter and Meadows ( 1995 1 to cover as large a range in (J - K) 

' IRAF is distributed by the National Optical Astronomical 
Observatories, which is operated by AURA, Inc., under cooperative 
agreement with the NSF. 



as possible, i.e. from -0.02 to 0.96, and they were observed 
at airmasses ranging from 1.03 to 1.80. Using the conversion 
from the Carter system to the CIT system determined by Carter 
( 1,1993 ) we find: 

^ciT = K,, + 0.014(+0.02) X (7 - r),„,, (1) 
-0.17(+0.03)xXh- const 

with an RMS scatter of 0.024 mag. We thus find no significant 
color term, and will assume a zero effect in the following. 

On the night of the IC4499 observations we also observed 
the standard star HD 136879 (Carter and Meadows T995^ be- 
fore and after the observations of IC4499 in a pattern of four 
positions on the detector array. The airmass for the standard 
star measurements was 1.44 and 2.13 respectively, while the 
airmass for the IC4499 observations ranged from 1.85 to 2.05 
(IC4499 is a very southern object). The standard star measure- 
ments were interpolated in airmass to provide the appropri- 
ate offset for the IC4499 data. Curves of growth were deter- 
mined from artificial aperture photometry of the standard star 
and an aperture radius of 20 pixels (10 arcsec) was used to de- 
termine the offset to the standard system. From the RMS scatter 
for Eq^we estimate the eiTor of the ^T-band zero-point to be 
0.025 mag. 

3.3. Photometry 

The DAOPHOT package (Stetson[T983 within IRAF was used 
to find the stars within the averaged frames and for performing 
photometry on them by point spread function (PSF) fitting. Due 
to the sparse nature of IC4499 crowding is not an issue for PSF 
fitting photometry. 

The aperture corrections for the PSF photometry were 
determined from a number of well exposed and reasonably 
isolated stars in the IC4499 fields. The aperture photome- 
try was performed after all the stars, apart from the stars se- 
lected for offset determination, had been subtracted by the 
IRAF/DAOPHOT substar routine. The resulting aperture 
corrections for five stars in each field agreed to better than 
0.02 mag. 

The RR Lyrae stars were identified from the finding chart of 
WN96. The final magnitudes on the CIT system are tabulated 
in Tabled together with the periods and the type of pulsation 
as determined by WN96. The pulsation type is indicated by 
"ab" for fundamental mode pulsators, "c" for first overtone, 
"d" for double mode, and "e?" for possible second overtone 
pulsators. The estimated errors provided by DAOPHOT have 
been tabulated as well. Note that the eiTors on the magnitudes 
for stars in field B are significantly larger than in field A due to 
the shorter integration time. 

The observed magnitudes have been collected at a specific 
time as already mentioned in SectlJ] so with an appropriate 
ephemeris for each star we can determine the phase of observa- 
tion. Alistair Walker has kindly provided the full optical light 
curves from WN96 in machine readable format and using their 
period estimates we have determined the epoch of maximum V- 
light for each of the stars. On this basis we have determined the 
phase corresponding to the observed /T-magnitude. We have 
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Table 1. Periods observed K magnitudes and the estimated 
mean K magnitudes for the RR Lyrae stars in the two fields. 
The fundamental period is listed for the RRd stars. The type of 
pulsator as determined by WN96 is listed as well. The last col- 
umn indicates in which of the two fields the star was observed. 



ID logP K a{K) (K) Type Field 

A~" 
B 
B 
A 
A 
B 
B 
B 
B 
A 
B 
A 
A 
A 
B 
A 
B 
A 
A 
B 
B 
B 
A 
A 
B 
B 
A 
B 
B 
B 
B 
A 

mag mag mag 



then used the ^T-band template light curves from Jones, Carney 
and Fulbright J1996> to compute the corresponding intensity 
averaged mean A'-magnitude, (K) , based on the type of pul- 
sation and the V-band amplitude as given by WN96. For the 
RRd stars we have not attempted a phase correction and for 
the putative RRe stars we have assumed that they are actually 
RRc stars as done by Castellani, Caputo and Castellani (2003 1. 
These estimates of the mean /T-magnitude have also been tab- 
ulated in Tabled 



4. Results 

4.1. The (Mk) -logP relation 

The (Mjf ) - log P relation has been established by Longmore et 
al. ( I1990> who determined the slope from photometry of RR 
Lyrae stars in 8 different galactic globular clusters. Recently, 



BOS has developed the relation further to incorporate a small 
but significant metallicity eff'ect and they find the relation: 

(Mk) = -0.770 - 2.101 logP + 0.231[Fe/H] (2) 

For a globular cluster like IC4499 where we assume that 
there is no significant spread in [Fe/H] among the RR Lyrae 
stars, we can simply insert the best estimate of the metallic- 
ity in Eq|2to get a (M/i-)-logP relation. WN96 discusses the 
available metallicity estimates for IC4499 and find a value of 
[Fe/H] = -1.65 + 0.1 based on spectroscopic measurements 
of Ca II triplet lines in four RGB stars by Cannon U992j and 
from a (V, (V - I)) color-magnitude diagram. 

We insert this value in Eq|2]and for each star we can then 
determine the distance modulus as: 

(m - M)k = {K) - {Mk) = W - 1.151 - 2.101 logPo (3) 

The relation refers to the fundamental period and thus it is 
necessary to fundamentalize the periods of the higher order 
pulsators. We adopt the period ratios as used by WN96, i.e. 
(P2:P1:P0)=(1.575:1.3439:1.0) and thus logP,, = logfi + 
0.128 and log Pq = log Pi + 0.197 for the RRc and RRe stars, 
respectively. 

In Fig[2we have plotted the {K) magnitudes for the indi- 
vidual stars versus log Pq, the fundamentalized period. Again 
we have assumed that the putative RRe stars are actually 
RRc stars. They actually fall exactly on the general relation. If 
we choose to fundamentalize the period of these stars as second 
overtone stars they would appear where the open diamonds in 
FigEhave been plotted which gives a less satisfying match. 

Inspecting Fig[2 further, we see that three of the stars, 
marked with open symbols (V14, V32 and V81), are overly 
bright for their period when compared to the bulk of the stars. 
WN96 also found both V32 and V8 1 to be over-luminous in 
the optical, and they suspected both stars to have an unresolved 
red companion. 

The detection limit is fainter than the faintest RR Lyrae star 
by about 1 mag in the shallow field (B) so the sample is com- 
plete and can not suffer from Malmquist bias. 

4.2. The reddening and the distance modulus 

For each star we can now use Eq|3lto determine the distance 
modulus and we can then compute the weighted mean of the 
resulting distances. If we use the sample of 29 RRab , RRc , 
and RRd stars excluding the three outliers we find (m - M)k = 
16.555 + 0.016 with cr = 0.087 mag. If we include V14 and 
V32 we find (m - M)k = 16.526 ± 0.023 with cr = 0.12 mag, 
which is only slightly smaller than the previous result. The 
three possible RRe stars do not affect the distance estimate sig- 
nificantly, independently of how we choose to fundamentalize 
them. Finally, if we choose to simply use the random phase K 
magnitudes as a proxy for the mean /T-magnitude instead of 
the template approach described above, we find (m - M)k = 
16.553 + 0.019 with cr - 0.104, basically an identical result but 
with a slightly larger uncertainty. 

Due to the relative insensitivity of the (Ma:) -log f rela- 
tion to reddening, we can determine an independent value for 



1 


-0.2147 


15.92 


0.06 


15.82 


ab 


4 


-0.2051 


15.79 


0.10 


15.80 


ab 


5 


-0.2542 


15.85 


0.10 


15.89 


ab 


11 


-0.1996 


15.73 


0.04 


15.79 


ab 


13 


-0.2911 


16.00 


0.05 


16.08 


ab 


14 


-0.3009 


15.67 


0.13 


15.63 


ab 


15 


-0.2242 


15.99 


0.12 


16.04 


ab 


17 


-0.3006 


16.05 


0.11 


15.96 


ab 


18 


-0.3248 


16.22 


0.14 


16.22 


d 


23 


-0.2944 


16.16 


0.07 


16.08 


ab 


24 


-0.2869 


15.98 


0.10 


16.06 


ab 


29 


-0.4405 


16.25 


0.08 


16.24 


c 


30 


-0.2766 


15.97 


0.05 


15.99 


ab 


32 


-0.4437 


15.33 


0.03 


15.31 


c 


37 


-0.2481 


15.82 


0.10 


15.87 


ab 


48 


-0.2842 


16.23 


0.08 


16.05 


ab 


49 


-0.3018 


15.91 


0.12 


15.98 


ab 


50 


-0.2474 


15.84 


0.06 


15.92 


ab 


51 


-0.3214 


16.20 


0.07 


16.20 


d 


55 


-0.4455 


16.08 


0.14 


16.10 


c 


56 


-0.5383 


16.21 


0.14 


16.24 


e? 


58 


-0.3005 


15.92 


0.12 


15.88 


ab 


59 


-0.3188 


15.97 


0.07 


15.97 


d 


70 


-0.2700 


15.98 


0.06 


15.87 


ab 


71 


-0.3045 


15.96 


0.11 


15.96 


d 


72 


-0.1666 


15.78 


0.09 


15.80 


ab 


81 


-0.4022 


15.71 


0.05 


15.70 


c 


82 


-0.2654 


15.83 


0.10 


15.88 


ab 


89 


-0.4493 


16.14 


0.14 


16.14 


c 


90 


-0.3113 


16.07 


0.13 


16.07 


d 


97 


-0.5402 


16.14 


0.14 


16.16 


e? 


1575 


-0.5534 


16.44 


0.09 


16.41 


e? 



4 



Jesper Storm: The distance to IC4499 



T 1 1 1 1 1 1 1 r 



15.5 



A 
V 



16 



16.5 




■0.4 



0.2 



-0.3 
log(Po) 

Fig. 1. The (K) magnitudes plotted against fundamental period. For the RRd stars the single epoch K is plotted. The open 
symbols were not included in the fit. The filled diamonds labeled RRe? represents the putative RRe stars from WN96 but 
fundamentalized as RRc stars. The open diamonds show where these stars would appear had they been fundamentalized as second 
overtone pulsators. The line shows EqlJlfor [Fe/H] = -1.65 and shifted by our best distance estimate of {m - M)k - 16.55. 



the reddening towards IC4499 from the mean V-magnitude of 
the RR Lyrae stars as determined by WN96 and the {My) - 
[Fe/H] relation. For consistency we use the calibrations pre- 
sented by BOS for both relations, and thus in addition to Eq. 
|2]we adopt 

(My) =0.177(+0.07)[Fe/H] + 0.718(+0.07) (4) 

which is vaUd for metal-poor stars ([Fe/H] < -1.6). 

WN96 find the mean V-magnitudes of the RR Lyrae stars to 
be< Vrr 17.652+0.006 and thus (m-M)y = 17.226+0.07. 
Adopting Ay = 3.1 X E{B - V) and Aj^ = 0.1 1 x Ay (CardelH 
et al. 1198 91 and using the constraint that the reddening cor- 
rected distance estimates from both the {Mk) -log P and (My) - 
[Fe/H] relations should be the same, we can write the relation: 



{m - M)v -Av-(m- M)k - Ak (5) 

from which follows: 

Ay = [(m-M)v-{m-M)K]/0.S9 (6) 

Ay = 0.755 + 0.08 (7) 

E(B - y) = 0.24 + 0.03 (8) 



Our best estimate of the reddening is then E(B - V) - 
0.24 + 0.03 and the best estimate of the distance to IC4499 is 
(m-M)o = 16.555-0.11x0.755= 16.472 ±0.016 where the 



error estimate is the formal fitting error only. The DIRBE red- 
dening map (Schlegel et alTgWl gives E(B - V) = 0.23 at the 
position of IC4499 suggesting that the reddening is completely 
dominated by galactic reddening. 

In addition to the fitting error we add the error contri- 
butions from the /T-band zero-point (0.025 mag), from the 
aperture correction (0.02 mag), from the uncertainty in Ak 
(0.009 mag), and from the uncertainty in the metallicity esti- 
mate (0.231 X 0.1 = 0.023 mag). This gives a final random 
error of 0.04 mag. The zero-point of the PL relation in Eq|2 
has an uncertainty of 0.044 mag which should be considered 
as a systematic error We should also keep in mind that the 
metallicity scale is not uniquely defined and that there might 
be systematic differences between scales reaching 0.2 dex 
(Layden '1994). This would constitute a systematic error of 
0.231 X 0.2 = 0.046 mag. This adds up to an overall system- 
atic error estimate of 0.06 mag. A much larger potential sys- 
tematic error arises from the adopted zero-point of the (Mk) - 
log P - [Fe/H] relation. As an example, Jones et al. (1992j found 
a zero-point based on Baade-Wesselink analysis of field RR 
Lyrae stars which is about 0.3 mag fainter than the B03 rela- 
tion. A similar result was found by Layden et al. (1996 J from a 
statistical parallax study of field RR Lyrae stars. Walker (I1992> 
showed that adopting the Baade-Wesselink and statistical par- 
allax zero-point would cause the LMC distance estimate from 
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RR Lyrae stars to be about 0.3 mag shorter than the canonical 
Cepheid based distance of 18.5 mag. This shows that the BOS 
relation is in good agreement with the canonical Cepheid dis- 
tance scale but it disagrees with the short distance scale. B03 
also show that their relation is in excellent agreement with the 
HST parallax measurements of RR Lyr itself (Benedict et al. 
[200 2 1. Kovacs 2003 has redone the Baade-Wesselink analysis 
for a sample of galactic Cepheids similar to the one of Jones 
et al. ( Il992t using the same calibration as has been success- 
fully used for Cepheids and he now finds agreement with the 
canonical LMC distance. We cannot explain the difference be- 
tween the distance scales with the current data, but we simply 
choose to adopt the B03 zero-point to ensure that we are on the 
canonical Cepheid distance scale. 

On this basis we find (m - M)o - 16.47 + 0.04 (random) 
+0.06 (systematic), or d = 19.7 + 0.4 + 0.6 kpc, as our best 
distance estimate to IC4499, which is in excellent agreement 
with the value of (m - M)o = 16.45 + 0.05 found by WN96. 
Our reddening value of E{B-V) - 0.24 + 0.03 also agrees well 
with the value of E(B - V) = 0.22 + 0.02 from WN96. 

We stress that for stellar systems containing a reasonable 
number of RR Lyrae stars with known periods and metallic- 
ities, the {MK)-logP -[Fe/H] relation combined with single 
epoch /T-band data gives accurate and reddening insensitive 
distances for a very modest amount of observational effort as 
the dominating source of error is not the number of individual 
phase points. 
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